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ABSTRACT 

We analyze observations obtained with the Chandra X-ray Observatory of 
bright Compton thick active galactic nuclei (AGNs), those with column densi- 
ties in excess of 1.5 x lO^^cm"^ along the lines of sight. We therefore view the 
powerful central engines only indirectly, even at X-ray energies. Using high spa- 
tial resolution and considering only galaxies that do not contain circumnuclear 
starbursts, we reveal the variety of emission AGNs alone may produce. Approxi- 
mately 1% of the continuum's intrinsic flux is detected in reflection in each case. 
The only hard X-ray feature is the prominent Fe Ka fluorescence line, with equiv- 
alent width greater than 1 keV in all sources. The Fe line luminosity provides 
the best X-ray indicator of the unseen intrinsic AGN luminosity. In detail, the 
morphologies of the extended soft X-ray emission and optical line emission are 
similar, and line emission dominates the soft X-ray spectra. Thus, we attribute 
the soft X-ray emission to material that the central engines photoionize. Because 
the resulting spectra are complex and do not reveal the AGNs directly, crude 
analysis techniques such as hardness ratios would mis-classify these galaxies as 
hosts of intrinsically weak, unabsorbed AGNs and would fail to identify the lu- 
minous, absorbed nuclei that are present. We demonstrate that a three-band 
X-ray diagnostic can correctly classify Compton thick AGNs, even when signif- 
icant soft X-ray line emission is present. The active nuclei produce most of the 
galaxies' total observed emission over a broad spectral range, and much of their 
light emerges at far-infrared wavelengths. Stellar contamination of the infrared 
emission can be severe, however, making long-wavelength data alone unreliable 
indicators of the buried AGN luminosity. 
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1. Introduction 

Accretion onto supermassive black holes powers active galactic nuclei (AGNs), but in- 
tervening material prevents direct views of most sources. The Compton thick AGNs, which 
have column density Nh > 1.5 x lO^^cm"^ along the hne of sight to the nucleus, can be 
particularly elusive, but these AGNs are astrophysically important. First, they are common, 
representing up to half of all Seyfert 2 galaxies (Risaliti, Maiolino, & Salvati 1999). Second, 
in synthesis models (e.g., Setti & Woltjer 1989; Comastri et al. 1995), they are essential to 
replicate the observed spectrum of the cosmic X-ray background, which peaks around 30 
keV (Marshall et al. 1980). 

The intrinsic X-ray continuum of an AGN up to about 100 keV may be described by 
a power law, typically of photon index F = 1.9, or spectral index a = 0.9 (Nandra & 
Pounds 1994). Intervening material preferentially absorbs soft X-ray photons, and the direct 
emission of the central engine is not detectable below 10 keV in Compton thick cases. In the 
obscuring medium, however, the intrinsic emission is strongly reprocessed. It emerges with 
an effective flat photon index and diminished intensity, less than about 1% of the intrinsic 
source luminosity. Thus, Compton thick AGNs easily fall below detection hmits of current 
X-ray surveys, although the intrinsic power of their central engines is not exceptionally feeble. 
These uncounted sources are responsible in part for the diminishing fraction of the X-ray 
background that is resolved toward higher energies (Worsley et al. 2005). 

In a Compton thick AGN, the large equivalent width (EW) of the Fe Ka fluorescence 
line is the most direct signature of its nuclear activity below 10 keV (GhiseUini, Haardt, 
& Matt 1994; Krolik, Madau, & Zycki 1994). The energy of this line is 6.4 kcV when the 
fluorescing medium is not fully ionized (less than Fe XVIII), or around 6.7 and 7.0 keV 
in Fe XXV and Fe XXVI, respectively. This line can also be an effective diagnostic of the 
geometry of the obscuring medium and the AGN's luminosity (Levenson et al. 2002). 

With the powerful continuum suppressed at lower energies, obscured AGNs reveal other 
soft X-ray spectral components. The active nucleus may produce line emission through pho- 
toionization and photoexcitation (Sako et al. 2000; Kinkhabwala et al. 2002; Ogle et al. 2003). 
The AGN may also drive an outflow that generates shocks in the surrounding medium, pro- 
ducing thermal emission. Alternatively, separate additional emission sources may be present. 
Speciflcally, compact circumnuclear starbursts are common in active galaxies (Cid Fernan- 
des, Storchi-Bergmann, & Schmitt 1998; Gonzalez Delgado, Heckman, & Leitherer 2001). 
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and they characteristically exhibit soft thermal X-ray spectra (Levenson et al. 2001), similar 
to starburst galaxies. Here we analyze seven high spatial resolution X-ray observations of 
optically-bright, Compton thick Seyfert 2 galaxies to identify X-ray and multiwavelength 
signatures of their obscuration and to measure the power of their hidden central engines. 

We define the sample in §2. We present the data and describe new results from the im- 
ages in §3. With the spectroscopy results (§4), we use the Fe K lines to determine properties 
of the buried AGNs, identify the soft X-ray emission as a consequence of photoionization, 
relate the soft X-rays to the continuum emission, and develop coarse X-ray discriminants of 
Compton thick AGN (§5). We discuss the reflection geometry of the continuum in §6. We 
use infrared (IR) and optical data to compare the observable and intrinsic characteristics of 
these galaxies across the electromagnetic spectrum (§7), and we summarize our conclusions 
in §8. 

2. Sample Members 

All sample galaxies are classified as normal Seyfert 2s on the basis of optical emission line 
ratios. We require that the sample members be previously identified as Compton thick and 
have bright AGNs, based on the core [O III]A5007 fine flux. The latter criterion allows us to 
measure the nuclear properties well. In order to reveal the variety of X-ray emission for which 
an AGN alone may be responsible, we exclude galaxies that are known to contain nuclear 
starbursts. Furthermore, to separate the immediate nuclear and circumnuclear environments, 
we require high spatial resolution, which only Chandra can provide. The exclusion of galaxies 
that contain starbursts and the spatial resolution requirement render this sample of Compton 
thick AGNs smaller than that of other recent work (such as Guainazzi et al. 2005a), but 
these restrictions allow us to isolate the active nuclei and therefore recover their intrinsic 
properties accurately. Data are then available for seven galaxies, down to core [O III] flux 
limit of 7 X 10~^^ergs~^ We rely on Whittle's (1992) [O III] fluxes, where available. We use 
the [O III] fluxes of ESQ 428-G14, NGC 3393, and NGC 5347 from Bergvall, Johansson, & 
Olofsson (1986), Diaz, Prieto, & Wamsteker (1988), and Gonzalez Delgado & Perez (1996), 
respectively. The sample is listed in Table 1. 

Koyama et al. (1989) measured a large EW Fe fine in a Ginga spectrum of NGC 1068 
and argued that the direct AGN continuum must therefore be completely blocked. Maiolino 
et al. (1998) identifled ESQ 428-G14, NGC 1386, and NGC 3393 as Compton thick in 
BeppoSAX observations. In the first, they noted the fiat (F = 0.57) spectrum, low ratio 
of observed 2-10 keV to [O III]A5007 fiux, and non-detection above 15 keV as indirect 
evidence for A^^;^ > 10^^ cm~^. In the second, they fit the 1-10 keV spectrum with a refiection 
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model and a prominent Fe Ka line. They fit the 1-80 keV spectrum of NGC 3393 to find 

Nh > 10^^cm~2. RisaUti et al. (1999) identified NGC 5347 as Compton thick on the basis 
of the large Fe Ka EW measured in a spectrum from ASCA and the low ratio of observed 
X-ray/ [O III] fiux. Risaliti et al. (2000) suggested that NGC 7212 is Compton thick from 
analysis of a low signal-to-noise ASCA spectrum, finding an apparently flat continuum and a 
prominent Fe line. Guainazzi, Matt, & Perola (2005b) confirm this identification in a higher 
quahty spectrum obtained with XMM-Newton. Awaki et al. (1991) first measured the large 
column density and prominent Fe K line of Mrk 3 in a spectrum they obtained using Ginga. 

Population synthesis analysis of optical spectra shows intcrmcdiatc-agcd and older stel- 
lar populations (ages > 100 Gyr) in the nuclei of nearly all sample members (Cid Fcrnandes 
et al. 2001, 2004, and references therein). The exception, NGC 5347, exhibits no significant 
emission from polycyclic aromatic hydrocarbons at 3.3//m (Imanishi 2003) and far-infrared 
(FIR) colors that are characteristic of AGNs (de Grijp, Lub, & Miley 1987). Thus, we do 
not expect any significant stellar contributions to the X-ray spectra of these galaxies. 

3. Observations and Image Analysis 

We observed the galaxies using the Chandra Advanced CCD Imaging Spectrometer 
(ACIS) back-illuminated S3 detector. We reprocessed all data from original Level 1 event 
files using Chandra Interactive Analysis of Observations (CIAO) software, version 3.2. (See 
the Chandra Science Center^ for details about Chandra data and standard processing pro- 
cedures.) We applied current calibrations (Calibration Database version 3.0), including 
corrections for charge-transfer inefficiency and time- dependent gain variations, to produce 
the Level 2 event files. In each case, we included only good events that do not lie on node 
boundaries, where discrimination of cosmic rays is difficult. We examined the lightcurves of 
background regions and excluded times of significant background flares. 

In the NGC 1068 observation, we use only the data with 0.1 s "frame times." With the 
CCDs read out after very short exposures, we minimize the problem of "pileup," in which 
coincident photons that arrive within a single CCD readout are not measured accurately. 
We otherwise applied the same data reduction procedures described above. This particular 
observation and observing mode required a correction to the on-source exposure^, for a net 
exposure of 1.53 ks. Mrk 3 was observed with the High- Energy Transmission Grating in place. 



^ http : / /cxc . harvard . edu/ 

^http: / /cxc. harvard.edu/ciao /caveats / acisjnterleave.html 
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We therefore do not analyze the extended emission of this galaxy, and we refer to published 
measurements of the spectrum (Sako et al. 2000) where comparisons are appropriate. 

Table 1 lists observation dates and resulting exposure times of the data we analyze 
directly along with some basic information about all of these galaxies. The Galactic col- 
umn density toward each source {Nh', column 4) is based on the Schlegel, Finkbeiner, 

& Davis (1998) measurements of Ay, with Nh = 1.9 x 10'^^ Ay cm~'^ mag~^ . We use 
Hq = 70 km s^^ Mpc~^ and redshift {z; column 5) to derive distances (column 6) and physical 
scales (column 7). 

The Chandra image of ESO 428-G14 (Figure 1) is typical. The images clearly show 
extended X-ray emission, demonstrating that the unresolvable central engine is not the sole 
observed X-ray source. The resolved emission extends over scales of hundreds of parsecs, 
and it is predominantly soft. We compare images at soft (0.3-1 keV), medium (1-4 keV), 
and hard (4-8 keV) energies with corresponding models of Chandra's point spread function. 
ESO 428-G14, NGC 1068, NGC 1386, and NGC 3393 are significantly extended in the soft 
and medium bands, even on small (1") scales, while their hard bands are unresolved. NGC 
7212 is significantly extended only in the medium band. The bright center of NGC 5347 is 
not resolved in any band, although this galaxy shows distinct emission approximately 3" to 
the northeast of the nucleus. 

In all cases, the large-scale extended X-ray emission appears similar to images of these 
galaxies in optical emission lines (Figure 2), as Bianchi, Guainazzi, & Chiaberge (2006) have 
shown for some of these sample members. For example, the morphology of ESO 428-G14 
in X-rays specifically and in detail resembles that of both Ha and [O III] (Falcke et al. 
1996). From the nucleus (evident at both radio and X-ray energies), very bright emission 
extends toward the northwest, with fainter emission extending toward the southeast. In 
Figure 2, we show the Chandra intensity contours overlaid on the continuum-subtracted Ha 
image obtained with the Hubble Space Telescope {HST) Planetary Camera of the Wide Field 
and Planetary Camera 2 (WFPC2). We aligned and scaled the image through the F814W 
filter (U2NP0603T) to subtract continuum from the F658N image (U2NP0601T). We then 
aligned the X-ray-identified nucleus (§4) with the maximum optical continuum emission. 
The maximum broad-band X-ray emission is located north of the Fe Ka emission peak that 
defines the X-ray center. Thus, the highest-intensity contours are not centered on the optical 
nucleus. Photoionization by the central engine produces the optical line emission in all these 
galaxies. While the lower spatial resolution X-ray observations cannot show variations on 
the very small scales (O'.'l) that are evident in the optical data, the similar overall X-ray 
morphology suggests that photoionization is also the origin of the soft X-ray line emission. 
The Chandra spectra also support this interpretation of the line emission as a consequence 
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of radiative processes rather than coUisional excitation as (§5.2). 

We find several sources in addition to the active nucleus near each galaxy's center. We 
present these results in Appendix A. 

4. Spectroscopy 

The ACIS data afTord spatially resolved spectroscopy of several interesting regions. We 
obtained spectra of the nuclei, extended circumnuclear regions, and the whole galaxies. We 
identified the locations of the AGNs in continuum-subtracted Fe Ka images. Line emission 
dominates the emission in images extracted with energy restricted 6 < £■ < 7 keV. We 
measure the corresponding continuum in the bandpass 4.5 < E < 5.5 keV, at which the 
ACIS response is similar. We smoothed each of these images by a Gaussian of FWHM = 1", 
then subtracted the continuum (4.5-5.5 keV) image from the corresponding linc+continuum 
(6-7 keV) image without any further scaling. This procedure is adequate to determine the 
locations of the maximum Fe Ka emission, but we do not use these images for quantitative 
measurements of the line flux. The hard X-ray peaks measured in the continuum-subtracted 
images are listed in Table 1. 

To isolate the AGN emission, we extracted spectra within an aperture of I'.'S radius cen- 
tered on each of these nuclear locations. Even in nearby galaxies, however, these unresolved 
"nuclear" apertures cover hundreds of parsecs around the active nuclei. While spectroscopy 
of large regions of some of these sample members has been presented elsewhere (Young, 
Wilson, & Shopbell 2001; Bianchi et al. 2006), the small nuclear apertures are critical in 
order to measure the immediate AGN emission, especially the Fe Ka line fiux and EW. In 
all cases, we subtracted the background measured in nearby areas that contain no obvious 
sources. Over the course of the Chandra mission, the soft X-ray sensitivity has diminished, 
likely the result of build-up of material on the detector. We created ancillary response files 
that account for this time- varying effect. We performed the model fitting in XSPEC (Arnaud 
1996). 

Broadly, the nuclear spectra are similar. The hard X-ray emission (E > 3 keV) is 
spectrally flat and relatively weak, with a prominent Fe Ka emission line in each case. The 
bright soft X-ray spectra show line emission. In order to measure the Fe K line accurately and 
to isolate the AGN continuum, we first fit the spectra from only 4-8 keV. We fit unbinned 
spectra using the C statistic (Cash 1979) to retain significant counts at energies greater 
than 6.5 keV. We model the Fe fine as an unresolved Gaussian, accounting for each galaxy's 
redshift in measuring the central energy of the line. 
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The observed flat hard X-ray spectra and previous evidence for Compton thick obscu- 
ration indicate that we cannot directly detect the intrinsic AGN power law continuum in 
these cases. Instead, we model the observed continuum as purely reflected AGN light in a 
neutral medium using the PEXRAV model (Magdziarz & Zdziarski 1995) in XSPEC. Within 
the Chandra bandpass, we are not sensitive to some of the model parameters (e.g., intrinsic 
photon index and high-energy cutoff), and others are degenerate. Thus, we fit only for nor- 
malization of the intrinsic power law, fixing the remainder of the model parameters: intrinsic 
photon index F = 1.9, high-energy cutoff = 300 keV, and abundances are solar (Anders & 
Grevesse 1989). The model describes reflection off a plane-parallel, semi-infinite slab, and it 
is a function of viewing angle, measured from the slab normal. We adopt inclination i — 63°, 
the model default. This is an intermediate value considering obscured views only; assuming 
a toroidal obscuring geometry with half-opening angle 6, the central engine is not viewed 
directly from i = 6 to 90°. Although the high-energy {E > 10 kcV) spectral shape varies 
with viewing angle, the primary effect on the Chandra spectra is to alter the normalization. 
Compared with the extreme values, the adopted inclination angle results in normalizations 
that may be overestimated by up to 25% (for true pole-on views) or underestimated by up 
to a factor of three (for equatorial views). We also account for Galactic absorption, although 
its effect is negligible. 

Table 2 contains the results of these hard nuclear fits, including the C statistic and 
number of degrees of freedom {dof). The tabulated values of the power law normalization 
(column 5) assume full (100%) reflection of the AGN continuum. With decreasing relative 
reflection, (e.g., with smaller covering fraction of reflecting material), the underlying AGN 
normalization increases proportionally. Figure 3 shows the hard spectral fits; the prominent 
Fe lines are evident in all these nuclear spectra. 

At the current sensitivity and resolution, there are no significant spectral differences 
among an intrinsic (heavily absorbed) continuum, this reflection model, and an unabsorbed 
F = power law. For example, the Fe absorption edge is strongest in the first of these 
models and absent in the last, but the edge cannot be distinguished in any of these data. 
The shortcoming of the flat power law is that this model is purely phenomenological; such 
intrinsic spectra are not observed in Seyfert 1 galaxies or less-obscured Seyfert 2s, where they 
would be detectable. The advantage of adopting the reflection model is that it does provide 
a lower limit on the intrinsic AGN power, which we will compare with other multiwavelength 
observations and calculations (§7). 

The soft X-ray spectra do not affect the high-energy results, so for simplicity, we fixed 
these latter model components while fitting the spectra over the full 0.4-8 keV range. We 
grouped the full spectra into bins having a minimum of 20 counts each and use the 
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statistic. While fixing foreground extinction at the Galactic values, we considered several 
physically-motivated model components to account for the soft X-rays: the MEKAL thermal 
model (Mewe, Gronenschild, & van den Oord 1985; Mewe, Lemen, & van den Oord 1986; 
Liedahl, Osterheld, & Goldstein 1995), additional power laws, and line emission. An AGN- 
driven outflow may generate thermal emission. The non-thermal continuum may be the 
intrinsic continuum reflected off the surrounding medium. Photoionization by the AGN 
continuum may produce recombination and fluorescence lines. 

Overall, wc find that the line emission dominates the nuclear spectra below 3 keV. 
While the ACIS detector cannot resolve individual transitions, the spectra show relatively 
sharp emission peaks, as opposed to the much broader blends of many lines that coUisional 
excitation would exhibit. We model these lines as unresolved Gaussians. In each spectrum, 
we include additional lines or other components as long as they are significant at the 95% 
confidence level, based on an F test. Figure 4 of ESO 428-G14 is representative of the 
full nuclear spectra and the best-fitting models. We emphasize that predominantly thermal 
models do not quantitatively fit any of the observed spectra as well as these line models do 
(§5.2). 

We consider the spectra of extended regions of the galaxies. (We exclude NGC 1068 from 
this discussion because its extended X-ray emission has been extensively analyzed previously; 
Young et al. 2001.) We extracted spectra from the bright, central regions of the galaxies, 
including their nuclei. These elliptical apertures are 6" x 3" in NGC 1386 and NGC 3393, 
and 5" x 3'.'5 and 3'.'5 x 2" in ESO 428-G14 and NGC 7212, respectively. Wc also obtained a 
spectrum of the soft X-ray emission northeast of the nucleus of NGC 5347, within a circular 
aperture of 1'.'5 radius. The total emission within 1' of each nucleus yields a spectrum of 
each galaxy as a whole. 

Tables 3 and 4 list the results of these model fits of the soft nuclear spectra and the larger 
central regions where our findings are new or significantly different from published values, 
and Table 5 summarizes the flux measurements. Some of the spectra do require model 
components in addition to the unresolved fines and hard X-ray features. The central ESO 
428-G14 spectrum requires an additional continuum source, which a power law describes 
well. Because the power law index is not well constrained, however, we flx its slope at 1.9, 
matching the intrinsic AGN continuum. 

Given the fimited spectral resolution of the detector, each modeled line in the soft 
spectrum represents a blend of emission lines and is not attributable to a single transition. 
However, we identify a particular element as the dominant emitter in several cases where the 
observed line energy corresponds to a known (measured or theoretical) value. For example, 
in all spectra we flnd emission near 0.56 and 0.91 keV, which we identify as O VII and Ne IX, 
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respectively. Several of the spectra suggest Si XIII near 1.8 keV. We do not characteristically 
find emission from fully-ionized species, although wc measure O VIII in NGC 1386 and it 
may appear blended in the spectrum of ESO 428-G14. Nc X may also contribute to the lines 
observed at 1.0 keV in several of the spectra. Most of the lines observed to have energies 
between 0.7 and 1.2 keV likely contain a strong contribution from the many Fe L shell lines 
in this energy range. 

The best-fitting models of the NGC 3393 nuclear and central spectra each include a 
thermal component. At this temperature, {kT ^ 1.6 keV), the collisional excitation results 
in a number of emission lines around 1 keV, where Fc L transitions arc also strong. We 
considered adding to the unresolved lines a broad Gaussian emission feature, to account 
for the many Fe L lines, instead of the thermal component. The model with the thermal 
component fits the data much better than this alternate purely descriptive model. In the 
nuclear spectrum, for example, x^/^ — 25/21 for Gaussian width a — 0.3 and central 
energy 0.9 keV. We note that both our hard and soft spectral fits differ from the results of 
Bianchi et al. (2006). The larger Fe Ka EW we find may be a consequence of our smaller 
aperture, selected to isolate the line emission against the AGN's continuum. We do not, 
however, significantly measure as many unresolved soft lines as they report, even in our 
larger aperture. 

The emission northeast of the NGC 5347 nucleus is distinct, but faint. Only the soft 
X-ray flux is significant. A single-temperature thermal model fits the unbinncd spectrum of 
this region, with kT = 0.24 ± 0.05, although optical line ratios indicate that this emission is 
photoionized (Gonzalez Delgado & Perez 1996). 

We fit the whole galaxy spectra with models based on the best-fitting nuclear and 
central models, freeing the both the hard and soft X-ray parameters. These spectra cover 
very large areas, so they contain the integrated flux of many different emission sources. We 
do not interpret the physical significance of the models in detail, but use these results for 
fiux comparisons (Table 5). 

Most of the total flux is located within the smaller (central and nuclear) regions, rather 
than extended throughout the host galaxies. The central regions account for more than 70% 
of each galaxy's total soft emission. Typically, the nuclear aperture alone contains most of 
the soft emission. NGC 3393 is exceptional; here less than 40% of the soft fiux is contained 
within the nucleus, and the larger (1700 x 770 pc) central area accounts for the plurahty of 
this galaxy's soft emission. In each case, the nucleus alone does provide at least half of the 
total observed galactic hard X-ray emission (measured within a 1' radius aperture). In most 
cases, however, the extra-nuclear fraction of observed hard X-rays is significant, around 40% 
to 50%. In a less-obscured active galaxy, nearly all of the hard X-rays would be detected in 
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the immediate vicinity of the unresolvable central engine, because the hard AGN continuum 
would be observed directly. In this group of Compton thick sources, the spectra show that 
the hard X-rays emerge indirectly. As a result, line emission at energies greater than 2 keV 
is significant and may arise on larger physical scales. Only in NGC 5347 is the net (line and 
continuum) reprocessing confined to a small region (radius < 300 pc). 

We estimate that X-ray binaries contribute a small fraction of the observed total hard 
X-ray flux of the whole galaxies we measured. Colbert et al. (2004) derive scaling relations 
for the 0.3-8 keV flux of unresolved sources in terms of galaxy luminosity measured at FIR, 
and UV wavelengths. To determine the UV luminosity, we use lUE observations of NGC 
3393 (Kinney et al. 1993). In the other cases, we very roughly estimate this value from 
observations at B, assuming f -If ° = 0.2 (Kinney et al. 1996). We find that the 0.3-8 
keV flux of X-ray binaries is less than about 10% of the total 2-10 keV flux. Modeling the 
typical source as a power-law (F = 1.8) spectrum with Galactic absorption, we conclude that 
the 2-10 keV flux of such point sources is less than 5% of the total 2-10 keV flux of each 
galaxy. 



5. Results 

5.1. Fe KcK Line Emission: An Indicator of the Intrinsic Nuclear Luminosity 

The line emission near 6.4 kcV is the most significant X-ray characteristic of the AGNs 
in these galaxies. All of the line central energies are consistent with Fe Ka in a "neutral" 
medium, i.e., less ionized than Fe XVIII. To remain neutral, the fiuorcscing medium must 
not receive the powerful ionizing flux of the central engine directly, and it is likely the torus 
of the unifled AGN scenario (Antonucci 1993), the same region that blocks broad optical 
emission lines from direct view in these Seyfert 2s. 

The line EW is a function of obscuring column density, covering fraction, and viewing 
angle (Krohk et al. 1994). Quahtatively, larger covering fraction produces larger EW, since 
more Fe-edge photons are captured and generate Ka photons. The continuum is most 
suppressed along the equatorial plane, so these viewing angles also result in larger EWs. 
The large equivalent widths conflrm the Gompton thick identification of these AGNs. At 
solar abundance, EWs larger than 1 keV require Nh > 1.5 x lO^^cm"^. 

The EWs place some limits on the covering fraction, dependent on the viewing angle, 
i (Krohk et al. 1994; Levenson et al. 2002). In these models, the intrinsic AGN spectrum 
has F = 1.9, and the obscuring torus has a square cross-section allowing for unobstructed 
views over half-opening angle 9. An EW < 1.5 keV does not strongly constrain the covering 
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fraction, yielding < 50°. With the large EWs of ESO 428-G14 and NGC 1386, we find 
^ < 35 and 20° for modest viewing angle (i = 65°), and 9 < A5 and 30°, respectively, viewed 
almost through the equatorial plane {i — 85°). 

The results of these numerical simulations also show that the Fe line luminosity and 
the intrinsic AGN luminosity are related. For each of the observed AGNs, we use these 
numerical models to find the region in inclination angle, covering fraction, and column density 
parameter space consistent with the measured EW and its errors. We then compute the mean 
ratio of Fe Ka to intrinsic continuum luminosity within this portion of parameter space. We 
typically find Lpe/ Lagn,2-w ~ 2 x 10~^, where Lagn.2-io is the AGN's intrinsic 2-10 keV 
luminosity. Table 6 lists Lpe (column 4) and the resulting values of Lagn,2-io (column 
5). We label this intrinsic 2-10 keV luminosity LAGN,Fe to distinguish it from the intrinsic 
luminosity we determine using other techniques. Typical of Seyfert galaxies, the outcome 
is LAGN,Fe ~ lO^^crgs"^. In all cases, the observed 2-10 keV luminosity represents only 
about 1% of the intrinsic luminosity {Lobs,2-io/ LAGN,Fe ~ 0.01.) Note that the observed 
luminosities Lobsfl.5-2 and -^065,2-10 in Table 6 are not corrected for any absorption; they 
correspond to the directly-detectable flux of the sources. 

Having the same underlying energy source, Lagn,2-w and L[oiii], the luminosity of the 
[O III]A5007 line, are correlated. Our Fe K estimates of the intrinsic luminosity arc consistent 
with this empirical relationship (Table 6, column 7). Specifically, we consider the Heckman 
ct al. (2005) result log(L4G'jv,2-io/-^[o///]) = 1-59 ± 0.48 for the optically-selected type 1 
AGNs. (These unobscurcd nuclei reveal the intrinsic relationship, without the significant 
correction for X-ray absorption that even the Compton thin type 2 AGNs require.) Figure 
5 illustrates the agreement of the Fe fine and [O III] estimates of intrinsic AGN luminosity. 
In addition to the sample members, which are uncontaminated by nuclear starbursts, we 
also plot the published results for two Compton thick galaxies that do contain starbursts 
(Levenson et al. 2004, 2005). The dotted fine is plotted at a luminosity ratio of one, not a 
fit to the data. 

We compare the NGC 3393 results with the BeppoS AX detection in the 15-220 keV band 
with the Phoswich Detector System (PDS) as a direct measurement of the intrinsic AGN. We 
use only the high-energy data to measure the intrinsic AGN continuum level, neglecting the 
lower-energy measurements, which include multiple emission sources. We use the pipeline 
processed data, fitting the standard grouped spectrum from 15 to 220 keV. Both power law 
and cutoff power law models fit the data, with F f» 2. Extending these models to the Chandra 
bandpass, we find Fagn,2-io = 6±2 x 10~^^ erg cm~^ s~^, or Lagn,2-io = 2±1 x lO^'^ergs"^, 
in reasonable agreement with the Fe line calculation. Maiohno et al. (1998) and Guainazzi et 
al. (2005a) also use the BeppoS AX data obtained with multiple instruments over a broader 
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energy range and obtain similar results. Most importantly, their analyses explicitly show 
the Compton thick obscuration and characteristic reflected spectrum that we detect below 
8 keV. 



5.2. Origin of the Soft X-rays: Extra-Nuclear Photoionized Gas 

The morphology (§3 and Bianchi et al. 2006) and spectra (§4) of the soft X-ray emission 
provide evidence that we observe photoionized material. The extended features of the X-ray 
images match those of optical images. In the spectra, we find that a limited number of 
strong lines produce the observed features, in contrast to broad line complexes that would 
be observed from thermal plasma at the current spectral resolution, as we demonstrate 
quantitatively here. Radiation (from the AGN) rather than collisions thus determines the 
ionization state of the emitting material. We therefore describe this gas as photoionized, 
although other radiative processes, including photoexcitation and resonant scattering, likely 
produce some of the observed lines. 

Figure 6 illustrates the general deficiency of the thermal models in describing these 
soft spectra, using ESO 428-G14 as an example. The best-fitting thermal soft X-ray model 
{kT — 0.70 keV) fits poorly {x^/u = 131/42). It results in a broad emission peak around 0.9 
keV, and it cannot reproduce the individual narrower peaks around 0.8, 0.9, and 1.1 keV that 
are evident at the present resolution. Two thermal components can fit the spectrum well; 
with kT = 0.63 and 2.1 keV, x^/'^ = 38/40. In this case, however, the difficulty would be 
to account for the hot component in a physically reasonable way. Shock velocities of nearly 
1500 km s~^ would be required, yet there is no evidence for the corresponding broad line 
widths in either these spectra or high-resolution spectra of other active galaxies, where they 
could certainly be detected. Purely thermal models produce similar results in the nuclear 
spectra of the other sample members. In summary, physically reasonable thermal models do 
not quantitatively fit the spectra well. 

Thermal emission can be identified and measured in CCD X-ray spectra when it is 
present. Starburst galaxies, for example, contain large volumes of hot gas, often in a galactic- 
scale wind (Heckman et al. 1990). Detailed analysis of Chandra ACIS observations of the 
starburst galaxy NGC 253 clearly shows thermal emission (Strickland et al. 2002). In this 
case, a two-temperature model (with kT — 0.2 and 0.7 keV) fits spectra of the diffuse 
emission well, although it requires low abundances. 

NGC 1068 serves as a useful example of photoionized material, in which radiative pro- 
cesses produce the soft X-ray fine emission. This brighter galaxy has been successfully 
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observed at high spectral resolution. Measurements of the strength of individual line tran- 
sitions reveals photoionization as the dominant excitation mechanism (Kinkhabwala et al. 
2002; Ogle et al. 2003). The corresponding spectrum obtained with ACIS offers a direct 
empirical model for the other spectra we discuss here. Figure 7 shows the spectrum from 
the central 1'.'5 around the nucleus of NGC 1068. It is similar to the others, showing rela- 
tively narrow features in the predominantly soft emission, which we fit well with multiple 
Gaussians. 

Wc measure more emission lines significantly in this high-quality spectrum compared 
with the the other sample members, and wc note that much of the flux at energies above 2 
keV is in lines, not continuum (Table 7). Young et al. (2001) analyze a spectrum extracted 
over a slightly larger (1'.'9) aperture in these same data. They show that the narrow line 
emission is significant, even when strong soft continuum components are included in the 
model. While interpretation of the physical nature of X-rays observed at low resolution does 
depend on the model-fitting procedure, we emphasize that high quality CCD spectra are 
sufficient to distinguish the thermal or radiative origin of the emission. 

Successful models of X-ray spectra of photoionized line emission exist, but we do not 
attempt to apply a more realistic model to these data for two reasons. First, although 
emission line ratios arc useful diagnostics of physical conditions, including identification of 
photoionized and collisionally-cxcitcd plasmas, we cannot measure any isolated lines in these 
data to distinguish among competing model descriptions. Second, the models arc sensitive 
primarily to the ionization parameter, the ratio of radiative flux to gas density, given the 
ionizing spectrum and distribution of illuminated material. In the centers of active galaxies, 
however, no single ionization parameter is appropriate. In such an inhomogeneous region, 
individual clumps of gas may have different densities or be located at different distances 
from the ionizing source, resulting in a range of ionization parameters. In fact, high-quality 
observations explicitly show a variety of ionization parameters within individual astronomical 
sources, with each line produced predominantly in the region in which its emissivity is a 
maximum (e.g., Sako et al. 1999). 



5.3. Relationship of Soft X-rays and Continuum Emission 

One important result of the spectral fitting is that the soft X-ray emission of the nu- 
clear and central regions consists almost entirely of lines. The soft X-ray emission of many 
Seyfert 2 galaxies has been attributed to the AGN continuum (e.g., Reichert et al. 1985; 
Awaki et al. 2000; Matsumoto et al. 2004). In this interpretation, although the central en- 
gine is generally obscured, the continuum could be observed through gaps in the covering 
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material. Alternatively, analogous to the broad emission lines that are sometimes observed 
in scattered (polarized) light, the obscured AGN continuum could be scattered into the line 
of sight. These partial covering or scattered continuum models for the soft X-ray emission 
were generally plausible in the past, given the capabilities of existing instruments. Viewed at 
low spectral resolution, lines appear blended into a pseudo-continuum, as others have noted 
while employing continuum models as phenomenological descriptions of their data (e.g.. 
Turner et al. 1997). However, the physical interpretation of these models — in terms of cov- 
ering fractions and scattering efficiencies, for example — will generate misleading conclusions 
if fundamentally lines comprise the modeled continuum. 

While all galaxies we analyze here exhibit strong soft X-ray emission, it is not the AGN 
continuum viewed in reflection or through patchy obscuration. Only the central spectrum of 
ESO 428-G14 includes such a scattered continuum component in the best-fitting model, and 
it accounts for less than one-third of the 0.5-2 keV flux. We consider the possible contribution 
of a scattered continuum component to the soft flux of all the nuclear spectra. Wc find strict 
upper limits of 10% in NGC 1068 and NGC 1386, 20% in NGC 3393 and NGC 7212, and 
26% in the nucleus of ESO 428-G14. Only the lower-quality spectrum of NGC 5347 does 
not provide a strong limit, allowing up to 60% of the soft flux to be scattered continuum if 
all the emission lines are much weaker than we find in best-fitting model. High-resolution 
spectroscopy of comparable AGNs yields similar results. For example, Kinkhabwala et al. 
(2002) detect no electron-scattered continuum in NGC 1068. The observed soft continuum 
of the Circinus galaxy accounts for less than 30% of the total 0.5-2 keV flux that includes 
well- measured fines in addition to the continuum (Sambruna et al. 2001). In Mrk 3, the 
continuum contributes less than 25% to the observed soft Chandra spectrum (Sako et al. 
2000), although Bianchi et al. (2005) indicate that the continuum accounts for over 50% of 
this galaxy's flux detected with the Reflection Grating Spectrometer of XMM-Newton up to 
1.5 keV, where the line emission is measured well. 

In general, the scattered continuum contribution relative to the intrinsic emission in- 
creases with decreased clumping factor and increased ionization parameter (Krolik & Kriss 
1995). As these authors demonstrate, scattering can also substantially enhance line emission. 
This effect is strongest for the resonance transitions, where the line scattered fraction can be 
a few times the average (continuum) scattered fraction. Resonance lines are most evident 
in the observed spectra, so scattering as well as direct recombination likely contributes to 
these features. 

We recognize that partial covering models may provide an accurate physical description 
of some other AGNs. The spectrum of NGC 4151, for example, has been interpreted in this 
context (Weaver et al. 1994; Ogle et al. 2000). Typically, the "obscured" sources in which the 



- 15 - 



nucleus appears partially covered at X-ray energies (e.g., Weaver et al. 1996; Immler et al. 
2003) exhibit at least some broad optical emission lines (classified as Seyfert types 1.5-1.9). 
The broad emission lines, which originate close to the central engine, are therefore not fully 
hidden from direct view. Provided that the covering fraction of the the broad emission line 
region itself is less than 1, partial unobscured views of the nuclear continuum are possible. 

Obscured active galaxies generally display soft X-ray emission greater than that of the 
absorbed AGN continuum. With their continua fully suppressed, these Compton thick AGNs 
show that the luminosity of photoionized line emission can be significant, and we suggest 
that it may contribute to such "soft excesses," in both obscured and unobscured galaxies. 
The 0.5-2 keV luminosity of the concentrated central regions is 0.3-11 x 10^°ergs~^, and the 
soft luminosity of the entire galaxies ranges from 0.4 to 14 x 10^°ergs~^. These values are 
comparable to the soft excess luminosity measured in other Seyfert galaxies (e.g., Turner et 
al. 1997), or the total soft X-ray luminosity of galaxies in which all the direct AGN emission 
is suppressed below 2 keV, i.e., where Nh > lO^^cm"^ (Guainazzi et al. 2005b). 

The soft X-ray emission of the Gompton thick AGNs we investigate here is not related 
in a simple way to their observed or intrinsic hard X-ray emission. The observed nuclear soft 
flux ranges from 0.1 to 1 times the observed hard X-ray flux, and the ratio of soft to hard flux 

in the central regions ranges from 0.1 to 0.8. Assuming the buried AGN generates the soft 
X-ray lines via photoionization, the ratio of observed soft to intrinsic hard emission indicates 
a reprocessing efficiency. This ratio varies from 1-3 x 10~^, relative to the Fe line estimate of 
Lagn,2-w- The production of soft X-rays through photoionization depends strongly on the 
properties of the surrounding material that may be illuminated, notably the optical depth 
and covering fraction. In NGG 5347, for example, the lack of dense material along direct 
lines of sight to the nucleus is likely responsible for the low reprocessing efficiency. 

5.4. A Diagnostic Diagram for Obscured AGN 

With limited data, such as in surveys, the observed ratio of soft and hard fluxes in 
the form of an X-ray hardness ratio is commonly employed to identify an AGN and to 
determine both the degree of obscuration and the intrinsic luminosity. The hardness ratio 
is an appropriate measurement if the buried AGN dominates the total X-ray emission and 
is detected directly in the hard X-ray bandpass. In such instances, the lack of soft X-ray 
emission accurately indicates the absorption while allowing direct measurement of the true 
AGN power at higher energies. In these Compton thick galaxies, however, the AGN is 
not observed directly, and the strong soft X-ray emission is not the AGN continuum. For 
example, defining the soft band flux {S) from 0.5 to 2 keV and the hard band flux [H) from 
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2 to 8 keV, the hardness ratio is {H — S)/{H + S). Evaluating the hardness ratio in the 

observed nuclear or galaxy spectra would incorrectly indicate that these are unabsorbed and 
intrinsically weak AGNs, having Nh < 2 x 10^^ cm~^ and intrinsic -^2-10 < 2 x lO^"*^ erg s^^. 
Figure 8 illustrates the agreement between the observed luminosity ratios and that of an 
unobscured AGN. 

The lack of moderate-energy (2-5 keV) photons in the presence of strong soft X-rays 
distinguishes the heavily obscured AGNs from weakly-absorbed yet genuinely faint sources. 
When the data do not allow spectral modeling, X-ray "colors" based on intensity ratios 
in at least three bands would help to identify these sources. Guainazzi et al. (2005b), for 
example, consider the effect of redshift and utilize a very soft (0.2-0.5 keV) band in their 
analysis. Alternatively, Levenson et al. (2005) demonstrate that large column densities can 
be identified when gas heated by stellar winds and supernovae produces the soft X-rays 
that contaminate the standard hardness ratio. We construct a similar grid of AGN models 
combined with with varying amounts of photoionization to explore the X-ray colors in soft 
(0.5-2 keV), medium (2-5 keV), and hard (5-8 keV) bands. We model the AGN continuum 
as either a F = 1.9 power law absorbed by 0-10^"^ cm^^ or the pure reflection continuum 
absorbed by 10^°cm~^ (for typical Galactic absorption). We separately consider including 
with the reflection model the flux of an Fe Ka line of EW= 1 keV. We adopt the empirical 
soft hne spectrum of ESO 428-G14 to represent the photoionized contribution. We mix the 
photoionized component with the various AGN models, scaling it to account for to 100% 
of the total observed counts with each continuum model. Finally, for all summed models, 
we compute the ratio of counts in the three bands Chandra would detect. 

Figure 9 shows the model intensity ratios as a "color-color" plot. The results can be 
ambiguous when the photoionized component accounts for a significant fraction of the total 
counts. In particular, three-band ratios alone do not distinguish between the pure reflec- 
tion model of the hard X-ray continuum and moderate absorption (Nh ~ 1-3 x 10^^ cm~^). 
With a significant fraction of photoionized counts in these Compton thin spectra, however, 
the photoionization efficiency is imrealistically large. For Nh = lO^^cm"^, a 50% pho- 
toionization contribution (by counts) requires that 6.5% of the intrinsic AGN's 2-10 keV 
luminosity emerge in lines in the 0.5-2 keV band. For 3 x lO^^cm"^, the corresponding 
luminosity fraction is 2.1%. In contrast, we observe a photoionized soft X-ray luminosity 
Lphotofi.b-2 < 0.005LAGAr,2-io in the high-quality spectra of the Compton thick galaxies of 
this study. We therefore suggest imposing an upper limit on the allowed photoionization 
efficiency. With the requirement Lphoto,'a.b-2/ Lagn, 2-10 < 0.01, the photoionized count frac- 
tion is restricted to be less than 14 and 33% for Nh — 1 and 3 x lO^^cm"^, respectively. 
Equivalently, this constraint rules out soft/medium count ratios above 0.2 and 0.9 in these 
two Compton thin cases. We conclusively identify Compton thick galaxies with a strict limit 



-17- 



of Lphoto / Lagn, 2-w > 0.03, the area to the right of the thick sohd hne in Figure 9. 
These excluded Compton thin mixed models are also marked in red. The model count ratios 
where 0.01 < Lphoto,o.5-2/ Lagn,2-io < 0.03 are plotted in green. 

Figure 9 also shows the observed Chandra count ratios of the nuclei of this study. The 
simplified approach does not reproduce the results of spectral modeling in detail, but all 
these AGNs would certainly be recognized as Compton thick. For example, NGC 1386 and 
NGC 3393 are somewhat weaker in the medium band than the template line emission, so 
their colors arc offset from those of the reflection model. In addition, we plot the observed 
colors of several contrasting Compton thin galaxies, selected from the available Chandra 
archive and not classified as the obscured counterparts of narrow-line Seyfert 1 galaxies, 
which may have unusual soft X-ray properties (Leighly 1999). The emission from NGC 5728 
was measured within a 2" radius aperture. The bright nuclei of NGC 2110, NGC 5252, and 
NGC 7172 are strongly piled up, so we calculated the count ratios in the resulting readout 
streak, which remains spectrally accurate. In the three-color diagram, all these contrasting 
galaxies would be classified as Compton thin, although none of their count ratios lie close 
to the values of the single-component absorbed power law model (at the extreme left of the 
solid curves). Finally, we plot the colors of the comparison Compton thick galaxies that also 
contain nuclear starbursts. Because the X-ray emission of the starbursts is almost entirely 
soft, the active nuclei in these composite galaxies remain correctly classified as Compton 
thick. We conclude that the combination of the three X-ray spectral bands and a realistic 
constraint on the efficiency can identify Compton thick sources without complete spectral 
modeling. 



6. Continuum Reflection Geometry 

The relationship between the observable refiected AGN continuum and the intrinsic 
emission is a function of geometry. Figure 10 schematically illustrates two scenarios that 
result in refiection-dominated spectra, assuming toroidal distribution of obscuring material 
around the central engine. The intrinsic AGN fight (solid lines) may be refiected off the 
far side of the torus and then viewed without obscuration (dotted line). Alternatively, the 
(near-side) obscuring material may simultaneously serve as the reflecting medium, in which 
case the reflected spectrum (dashed line) is also absorbed, by an amount that depends on the 
depth within the material at which it is produced. In cither case, both the heavily-absorbed 
intrinsic emission transmitted through the obscurer and the reflected component contribute 
to the net observed spectrum, but the transmitted component is weak in the Compton thick 
sources. For simplicity, the cartoon shows a uniform obscurer, although more realistically. 
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the torus may be inhomogeneous. A clumpy torus results in the same X-ray spectra, provided 
that the total obscuration along the line of sight is unchanged and — in the first scenario — 
individual clumps remain optically thick to act as reflectors. (Indeed, rapid X-ray variability 
observed in other obscured AGNs, such as NGC 1365 (Risaliti et al. 2005), supports the 
suggestion of a non-uniform torus.) The cartoon also shows the soft X-ray emitting region, 
which has an unimpeded view of the AGN continuum. 

Unlike less-obscured examples, in which some portion of the AGN emission is observed 
directly, recovering the AGN luminosity from the modeled continuum of these Compton 
thick galaxies docs not yield certain results. Interpreting the model fits to the spectra in 
terms of either of these scenarios yields a lower limit on Lagn,2-io- In the context of the first 
scenario, we assume the ratio of reflected to intrinsic emission is 1, then compute the 2-10 
keV luminosity of the unseen (buried) continuum source. The resulting values of Lagn,2-io 
are listed in Table 6 (column 6, labeled Lagn,pex)- Because the relative reflection is in fact 
less than 1, with some covering fraction less than 100% by necessity allowing direct views of 
the reflected continuum, this method gives a lower limit on Lagn,2-iq- 

This estimate of Lagn,2-iq is also a lower limit if the second scenario is applicable. As 
in the first situation, the covering fraction of reflecting material may be less than 1. The 
reflected spectrum may further be obscured by the foreground layers of the torus. Viewed 
through a screen of Nh = 10^"^ cm^^, for example, the emerging flux at 8 kcV is reduced 
by a factor of about 5. Indeed, the values of Lagn,2-v) based on the continuum modeling 
are about an order of magnitude lower than the reasonable estimates wc derived from the 
Fe line luminosities (§5.1), with the exception of NGC 7212. Similar reflected and absorbed 
descriptions of the observed X-ray emission of obscured AGNs have been proposed previously. 
Awaki et al. (2000), for example, applied such corrections to a sample of Seyfert 2s and 
increased the estimated intrinsic hard X-ray luminosity by factors of up to 10 over a single 
foreground absorber model. 

The reflected soft X-ray emission is detectable in the flrst case, while it is absorbed 
in the second case. The photoelectric cutoff energy (below which the continuum does not 
emerge) is a function of the overlying column density, and a cutoff at high energy (indicating 
large column density) would arise only in the second scenario. While the resulting AGN 
spectra differ, we cannot discriminate between these situations in the data. Line emission 
obviously dominates the very soft X-ray spectra, so the presence or absence of the soft 
reflected continuum is uncertain. At higher energy, where the spectral cutoff could decidedly 
identify the second scenario, the limited resolution and sensitivity prevent distinction of a 
true continuum from blended lines. 

As long as some portion of the far side of the torus is visible without obscuration. 
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however, this reflected light dominates the observed hard emission. This flrst scenario is 
more important even when the viewing angle is significantly greater than the torus opening 
half-angle. For example, with a viewing angle i = 70° and 6 = 40°, the 8 keV flux of 
the far-side reflection is nearly an order of magnitude greater than the emerging near-side 
reflection. In the present data, the existence of the photoionized regions demonstrates that 
the covering fractions of the tori are less than 1, so the flrst scenario best describes the origin 
of the observed reflection spectra. 

7. Multiwavelength Comparisons 

While deeply buried, the active nucleus remains the primary energy source in each of 
these galaxies. We compare the X-ray results with several other indicators of the luminosity 
of these very obscured AGNs at other wavelengths. Much of the intrinsic luminosity emerges 
at IR wavelengths, but stellar emission is a source of confusion. The extreme ultraviolet 
(EUV) dominates the intrinsic spectral energy distribution (SED), and we probe this regime 
indirectly through optical recombination lines. We demonstrate that the intrinsic properties 
of these AGNs are similar to their less obscured and unobscured counterparts. However, 
their observable characteristics are distinct, and therefore offer some identifying signatures. 

These Compton thick galaxies generally exhibit 25- to 60-//m flux ratios typical of 
AGNs (/25//60 > 0.26; de Grijp et al. 1987), in contrast to cooler, stellar- heated dust. We 
calculate the total 8-1000/im infrared luminosity from IRASHuxes, following the prescription 
of Sanders & Mirabel (1996), and flnd Lm ~ 10"^^ergs~^ (Table 6). The total IR luminosity 
indicates the bolometric luminosity; in unobscured quasars, Lhoi ~ 5Ljf>. Ignoring any 
stellar contributions and applying this relationship to these Compton thick sources suggests 
Lboi ~ 10^^-10^5 erg here. 

While the ratio of directly-observed X-ray to Lj^ is very low for AGNs, with Lobs,2-io/ Lir < 
10~^ in all cases, considering the intrinsic X-ray luminosity, LAGN,Fe/ Lir ~ 0.1. These latter 
flux ratios arc similar to those measured in unobscured sources. For example, in the median 
radio-quiet quasar spectrum of Elvis et al. (1994), Lagn,2-w/ Lir = 0.2. Risaliti & Elvis 
(2004) derive the average quasar SED in several broad bands, including corrections for sam- 
ple selection effects. This average spectrum has Lagn,2-io/ Ljr — 0.05-0.2, depending some- 
what on the assumed hard X-ray spectral shape and treatment of the shortest-wavelength 
IR emission. 

Neither the composite SEDs nor our AGNs have been corrected for star-heated dust 
emission at IR wavelengths, which can be important within the large IRAS apertures. For 
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example, while the kpc-scale star-forming ring of NGC 1068 is easily excluded from the 
nuclear X-ray measurements, it does contribute significantly to the IR measurements. The 
comparison of this sample with active galaxies that also exhibit significant star formation 
demonstrates that L/^ is not a reliable indicator of intrinsic AGN luminosity when the stellar 
populations are unknown. While the correlations of IR luminosity with AGN indicators 
based on [O III] and Fe line luminosity are similar for the non-starburst Compton thick 
AGNs and unobscured AGNs, the active galaxies that also contain strong star formation 
show IR luminosities that are up to an order of magnitude higher (Figure 11). 

Mid-infrared emission requires hotter dust, so the stellar contribution to the integrated 
luminosity is proportionally smaller in this restricted wavelength range (8-40/xm). As Figure 
12 illustrates, however, even the mid-IR excess of the star-forming galaxies is significant. 
Similar to Figure 11, the non-starburst Gompton thick galaxies exhibit luminosity ratios 
comparable to those of unobscured AGNs. Although the mid-IR does not always dominate 
the total IR luminosity, the two are very well- correlated in non-starburst galaxies because 
the underlying energy source — the AGN — is the same. 

The large-scale optical emission provides one more method for constraining the intrinsic 
luminosity of these AGNs. The recombination rate of hydrogen is related to the Ha hne 
luminosity A^^ec — 7.4 x IO^^Lho, where Lna is measured in erg s~^ (Mulchacy ct al. 1994). 
Assuming every hydrogen- ionizing photon that encounters circumnuclear material eventually 
results in recombination, the ratio of the recombination rate to the total ionization rate is 
the covering fraction, C, of photoionized material. We assume a power-law spectrum from 
the optical to X-ray bandpass, Li^ oc z/^", with a = 1.5. Integrating this continuum from 
the hydrogen ionization edge through 2 keV, the ionizing luminosity Lion oc LuaC'^. 

We measure Lu^ and C from the continuum-subtracted HST image of ESO 428-G14. 
In the resulting image, we find that line emission extends over the 3'.'5 x I'.'l region around 
the nucleus, for C = 0.20. We use the WFPC2 exposure time calculator'' to convert the 
observed count rate to flux for the line spectrum. We correct for the [N II] emission which 
accounts for about 1/3 of the flux in the fllter bandpass here (Falcke et al. 1996). We also 
correct for the Galactic extinction to yield Fh^ = 8.1 x 10~^^ergcm~^s~^. 

Ferruit, Wilson, & Mulchaey (2000) report a flux of 5.4 x 10~^^ erg cm"^ s"-*^ from Ha 
and [N II] in a 6" x 8" aperture at the center of NGG 1386. Because this measurement 
was made through HSTs hnear ramp fllter, the sensitivity to the [N II] doublet is less 
than that of Ha. We do not correct for the [N II] contamination here. We measure the 
covering fraction of the extended Ha emission in the PG image, flnding C = 0.17. Gooke 



''http://www.stsci.edu/instrumeiits/wfpc2/Wfpc2_etc/wfpc2-etc.html 
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et al. (2000) report a flux of 3.67 x 10"^^ erg cm-^ g-^ from Ha and [N II] in a 7'/5 x 12'/4 

area around the nucleus of NGC 3393. Correcting for Galactic extinction and an estimated 
33% contamination of [N II], Ff/Q, = 4.4 x 10~^^ erg cm~^ s~^. We measure the covering 
fraction of the S-shape evident in the i/^T continuum image, which line emission dominates, 
and find C = 0.26. We determine the covering fraction of NGC 5347 in the continuum- 
subtracted [O III] image from HST (Schmitt et al. 2003). Measuring the opening angle to 
the northeastern emission, we find C = 0.22. In narrow-band images, Gonzalez Delgado & 
Perez (1996) find Fffa — 7.1 x 10~^^ergcm~^ s~^ from the nucleus and northeastern knot 
together, after correcting for [N II] emission and extinction. 

Table 6 lists L^a and the calculated values of I/ion! the latter are generally around 
lO^^ergs"^, except for the lower-luminosity NGC 1386. This method of calculating Li^n 
yields a lower limit on the true value. The large-scale geometry of the emission-line region 
determines the value of C. Because the gas is clumpy on small scales, however, the true cov- 
ering factor could be substantially smaller, resulting in larger values of Lj„„. Also, while we 
correct the line fiuxes for Galactic extinction, we do not correct for intrinsic dust absorption. 

We associate the X-ray photoionized material with the optical narrow line region. The 
observed line luminosity Lna is typically an order of magnitude larger than the total observed 
0.5-2 keV luminosity. The total luminosity of the narrow line region in optical and UV lines 
is then a few xlO^^ to lO'^^ergs"^, or several orders of magnitude greater, even without 
correcting for intrinsic extinction or covering fraction. Fundamentally, these luminosity 
difi^crcnccs demonstrate that the emission region is optically thin to X-rays but optically 
thick in the extreme ultraviolet (EUV), the ionizing continuum of the optical lines. 

Another way to quantify this difference is in terms of the emission measure, EM = 
J n'^dl, where n is density, integrated along the line of sight. Wc treat all of the emission 
in the X-ray "lines" near 0.56 and 0.91 kcV as O VII and Nc IX, respectively, using the 
maximum line emissivities from the XSTAR photoionization code (Kallman & Bautista 
2001). Table 8 fists the emission measures, evaluating over the central regions of ESO 428- 
G14, NGC 1386, and NGC 3393, and the nucleus of NGC 5347. We use the Ha areas and 
luminosities above to find EM Ha- The X-ray emission measures are generally several orders 
of magnitude smaller than the corresponding EMna, even though wc overestimate all the 
X-ray line fluxes by counting blends of many unresolved lines. The Ha emission measure of 
NGC 3393 is relatively low, a consequence of the very large area (7'.'5 x 12'.'4) over which the 
luminosity was measured, while the bright Ha emission is confined to a much smaller area 
(2" X 4"). 

The viewing geometry of these Seyfert galaxies reveals the X-ray line emission and 
optical narrow line regions directly, and these may be related to the "warm absorbers" and 
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UV absorbers observed in unobscured AGNs (Blustin et al. 2005; Crenshaw et al. 2003). The 

present observations and known warm absorbers are different in several ways, however. First, 
the X-ray emission measures indicate somewhat lower column densities (e.g., Nqvii ~ 10^* 
rather than lO^^cm"^), although these are uncertain by an order of magnitude, given the 
unknown covering fractions. Second, the lower ionization states, such as O VII and Ne IX, are 
the most prominent emission features here. (High resolution spectra reveal some emission 
from fully-ionized species, but it tends to be weaker and have smaller emission measure; 
Sako et al. 2000; Kinkhabwala et al. 2002.) In contrast, fully-ionized species predominate 
in the absorption spectra (e.g., Netzer et al. 2003; Yaqoob et al. 2003). Finally, the spatial 
scales of the emission and absorption appear to be different. While we measure emission 
on scales of hundreds of pc, UV and X-ray absorbers are located within tens of pc of the 
AGN (Crenshaw & Kraemer 2005, and references therein). Although they are not identical, 
the absorbing and emitting material need not have entirely distinct origins to account for 
these observed differences. Instead, each preferentially samples a common distribution of 
properties differently: "warm absorption" spectra reveal the smaller scale, rapidly outflowing 
material, while the material located farther from the nuclei dominates emission. 

The multiwavelength data overall answer two important questions: How is the bolomet- 
ric luminosity of the AGNs distributed across the electromagnetic spectrum, and at which 
energies are the active nuclei revealed? Most of the luminosity clearly emerges in the IR. 
The IR spectra and near-constancy of intrinsic X-ray to IR luminosity ratios demonstrate 
that the buried AGNs produce the observed IR emission in these galaxies lacking starbursts. 
With our calculation of Ljon, we indirectly recover most but not all of the intrinsic EUV 
luminosity. The AGNs themselves preferentially emit in the EUV, and this is the dominant 
bandpass of quasar spectral energy distributions (Elvis et al. 1994), although the emergent 
EUV radiation of all obscured AGNs, even the Compton thin variety, is weak. Indeed, we 
find rough agreement between and as estimates of the bolometric luminosity of 
each of the buried AGNs. The sense of the discrepancy between the two wavelength regimes 
{Lion < Lir) is expected, for the calculated Lion is a lower limit and the observed Lm may be 
contaminated by star-heated sources in the large aperture. We conclude that the underlying 
spectra of these Compton thick AGNs are similar to those of unobscured sources, although 
the observable spectra are very different. 

The AGN characteristics of Compton thick sources are not different from other obscured 
AGNs at optical wavelengths, so there is no special bias against finding them in surveys. 
The general disadvantages of optical identification of Seyfert 2s are relevant, particularly 
the necessity of obtaining spectra rather than images alone, and host galaxy dilution of the 
AGN spectral features, which becomes more severe at high redshift (Moran, Filippenko, & 
Chornock 2002). 
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X-ray surveys efficiently find active galaxies, despite the small fraction (< 5%) of the 
total luminosity even unobscured AGNs emit in this bandpass. The high sensitivity of the 
current detectors and the high contrast with stellar sources (which are proportionately even 
fainter in X-rays) particularly contribute to this success. The Compton thick examples here, 
however, illustrate the types of sources that X-ray surveys can miss. 

X-ray observations with even modest spectral resolution do offer one technique to iden- 
tify Compton thick AGNs specifically: the large EW of Fe Ka, as we and others have 
previously suggested (Levenson et al. 2002; Maccacaro et al. 2004). Measuring the line 
against the suppressed continuum docs not require complex spectral modeling. The line 
emission is confined to a narrow spectral bandpass, so the desired signal may easily stand 
out against the low background of Chandra and XMM-Newton within a limited energy range. 
Furthermore, with the EW as a constraint, the line luminosity is an effective indicator of the 
intrinsic AGN luminosity. 

Compton thick AGNs may reveal themselves in the IR, having long- wavelength spectral 
slopes similar to other types of AGNs. An important caveat for IR identification, however, 
is that the AGNs do not always dominate the IR emission, even at mid-IR wavelengths. 
Starbursts can be powerful, resulting in IR spectra that do not indicate the active nuclei 
that are present. Such stellar contamination is a real concern, because starbursts are common 
in active galaxies. 

The Compton thick AGNs themselves may be viewed directly only at high energy (E > 
10 keV), yet even then obscuration in excess of lO^^'^cm"^ can reduce the emergent high- 
energy fiux by two orders of magnitude. Provided the obscuration is not so severe, the 
SWIFT Burst Alert Telescope survey should detect Seyfert galaxies like these and provide 
direct measurements of the intrinsic AGN luminosities. The limiting 15-150 kcV flux of the 
survey will be 1 mCrab after two years of operation (Parsons et al. 2004). Using the Fe 
line-determined Lagn,2-io to predict the harder emission, these sources have flux around 
2-3 mCrab. Other bright Seyfert galaxies with intrinsic Lagn,2-iq — lO^^ergs"^ could 
be detected to a distance of nearly 80 Mpc. Truly luminous obscured AGN having X-ray 
luminosities two orders of magnitude greater would be detectable to nearly z — 0.2. 



8. Conclusions 

Nearby Seyfert galaxies serve as important examples of AGNs in general and Compton 
thick sources specifically. Because stellar sources do not significantly contaminate the X-ray 
emission and the direct AGN fight is blocked, the present observations sensitively reveal 
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the variety of X-rays that AGNs alone can produce. While the underlying continuum may 
be a simple power law, the radiation that emerges is reprocessed, in soft photoionized lines 
and a weakened hard continuum. We find no evidence here that the soft X-ray emission is 
primarily the AGN continuum viewed through patchy obscuration or scattered into the line 
of sight. When active nuclei are viewed directly, as in type 1 sources, the powerful continuum 
dominates the entire observed spectrum. With the direct soft X-rays blocked, however, even 
by smaller column densities than in these examples, the detectable soft X-rays may be in 
emission lines. 

These spectra clearly show that line emission, not a smooth continuum, composes the 
soft X-ray spectrum. Although we do not spectrally resolve individual transitions, the present 
resolution is sufficient to identify photoionization as opposed to coUisional excitation as the 
origin of the hne emission. In all cases, the correlation of the extended soft X-rays with 
photoionized line emission at optical wavelengths supports this interpretation. 

Because the observed soft X-ray emission is not the AGN continuum viewed directly, 
the commonly-employed hardness ratio does not accurately measure either the obscuration 
or the intrinsic luminosity of these sources. Instead, we propose a three- "color" diagnostic 
to characterize AGNs — including the Compton thick ones — when high-quality data are not 
available for spectroscopy. 

Typical of Compton thick AGNs, even the hard X-rays that Chandra detects here are 
reflected, being only about 1% of the intrinsic AGN power. The clearest X-ray signature of 
the active nucleus in each of these galaxies is the large EW Fe Ka line emission, with EW 
> 1 keV in all cases. The EWs are large because the hne is measured against a reflected 
continuum. 

We determine the intrinsic luminosity several ways. The most reliable method in these 
data comes from the Fe line luminosity. The hidden emission emerges in the IR, and the 
intrinsic X-ray/lR flux ratios are typical of unobscured quasars, except where stellar IR 
contributions are strong. Because star-heated dust may be luminous, even at mid-lR wave- 
lengths, IR data alone do not accurately measure the AGN luminosity. The X-ray continuum 
modeling provides only a lower limit on the AGN luminosity, independent of the assumed 
reflection geometry. The reflection fraction is unknown, and the reflector may itself be 
obscured. 

This research has made use of the NASA/IPAC Extragalactic Database (NED) which 
is operated by the Jet Propulsion Laboratory, Cahfornia Institute of Technology, under con- 
tract with the National Aeronautics and Space Administration, and it has made use of data 
obtained from the High Energy Astrophysics Science Archive Research Center (HEASARC), 
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provided by NASA's Goddard Space Flight Center. Some of the data presented in this pa- 
per were obtained from the Multimission Archive at the Space Telescope Science Institute 
(MAST). STScI is operated by the Association of Universities for Research in Astronomy, 
Inc., under NASA contract NAS5-26555. This work was supported by NASA through grant 
G04-5117 and NSF award AST-0237291. NAL thanks the GSFC X-ray Astrophysics Labo- 
ratory for their hospitality and the NSFRO for support during part of this work. 

Facilities: BeppoSAX (PDS), CXO (ACIS), HST (WFPC2) 

A. Additional Sources 

In the Chandra images, we identify a number of sources in addition to the AGNs within 
2' of each galaxy's center. The sources are detected with a minimum of 3 net counts and 
non-zero spatial extent using the GIAO wavdetect algorithm. We considered total (0.3-8 
keV), low-energy (0.3-2 keV) and high-energy (4-8 keV) images separately, but the sources 
of the latter two images are a subset of sources in the total-band image. In all cases, we 
measure the counts in a 3" radius aperture after subtracting a local background. The source 
locations and their net count rates are listed in Table 9. The table also contains the projected 
angular distance of each source from the nucleus of the Gompton thick galaxy in the field 
and the corresponding physical scale at the distance of the central galaxy. 

We very crudely estimate the observed fluxes of these sources, assuming that their 
spectra are power laws with F = 1.9, absorbed by the Galactic foreground column density. 
We find that 1 x 10~^ count corresponds to 8 x 10~^^ ergcm"^ s~^. Individual sources 
with total X-ray luminosity > 10^^ erg s~^ are of interest as "ultraluminous X-ray" candidates 
(Miller & Golbert 2004). This limit corresponds to observed count rates of 17, 67, 4, 10, 
and 0.8 xlO"^ counts 8"^ in ESQ 428-G14, NGC 1386, NGC 3393, NGG 5347, and NGG 
7212, respectively, assuming the sources are located in the host galaxy. However, the large 
apertures in which the source counts are measured here encompass hundreds of parsecs and 
therefore likely include many distinct sources. 

These additional sources may be entirely independent of the galaxies at the centers of 
these fields. The sources around ESQ 428-G14 and NGC 1386 tend to be soft, however, and 
the number of sources in each of these region exceeds the source count rates of blank sky 
fields. Applying the results of Tozzi et al. (2001), for example, we would expect fewer than 
about 5 sources with 0.5-2 keV count rates > 1 x 10~^ count s~^ in each of the 13n' areas we 
surveyed, while we find 17 and 13 around ESQ 428-G14 and NGC 1386, respectively These 
facts suggest that many of the sources are indeed associated with (and possibly located in) 
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these host galaxies. Because the total number of sources is small, there is no significant 
trend in the number of sources as a function of distance from the center of each galaxy. The 
number of sources we find near NGC 3393, NGC 5347, and NGC 7212 is not significantly 
greater than the detection rates expected from blank field observations. We do measure soft 
emission from all of the additional sources near NGC 3393 and NGC 7212, and from half 
the sources near NGC 5347. 
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Table 1. Observations and Galaxy Data 



Galaxy Observation Date Exposure (ks) Nh (cm ^) z Distance (Mpc) scale (pc/") RA'' Decl* 



ESQ 428-G14 


2003 Dec 26 


29.6 


1.2 X 10^1 


0.0057 


24.3 


120 


07 16 31.18 


-29 19 28.92 


NGC 1068 


2000 Feb 22 


1.53 


2.2 X 10^0 


0.0038 


16.2 


80 


02 42 40.69 


-00 00 48.11 


NGC 1386 


2003 Nov 19 


19.6 


7.7 X 10^^ 


0.0029 


12.4 


60 


03 36 46.18 


-35 59 57.33 


NGC 3393 


2004 Feb 28 


28.0 


4.7 X lO^f 


0.0125 


53.6 


260 


10 48 23.44 


-25 09 43.53 


NGC 5347 


2004 Jun 05 


33.6 


1.3 X 10^0 


0.0078 


33.4 


160 


13 53 17.80 


-1-33 29 27.41 


NGC 7212 


2003 Jul 22 


19.9 


4.5 X lO^o 


0.0266 


114 


550 


22 07 02.00 


+10 14 01.0' 


Mrk 3 






1.2 X 10^1 


0.0135 


57.9 


280 







Measured from peak Fe emission. 
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Table 2. Hard X-ray Spectral Parameters 



Galaxy pe pe EW^ Pe agn C/dof 



ESO 428-G14 


6.35 ±0.03 


7.6 ±2 


1.6 ±0.5 


1.7 ±0.3 


238/270 


NGC 1068^ 


6.42 ±0.05 




1.2 ±0.7 


15 ±0.5 


115/201 


NGC 1386 


6.39 ±0.03 


5.5t^ 


2.3 ±1.5 


0.8 ±0.2 


173/270 


NGC 3393 


6.41 ±0.07 


4.2 ±2 


1.4 ±0.7 




189/270 


NGC 5347 


6.38 ±0.03 


4.7 ±2 


1.3 ±0.5 


1.2 ± 0.2 


231/270 


NGC 7212 


6.39 ± 0.04 


8.1 ±3 


1.0 ±0.4 


2.7 ±0.5 


277/270 



Note. — Errors are 90% confidence limits for one interesting parameter. 
'^Energy of line center in keV. 
''Line flux in 10~^ photons cm~^ s~^. 
'^Equivalent width of line in keV. 

''Normalization of intrinsic power law in units of 
10~^ photons keV~^cm~2s~^ at 1 keV. 

'^Fit over 5-8 keV. 
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Table 3. ESO 428-G14 Soft Spectra 



Nucleus Center 

-'^ une une -'^ une ^ line 



0.61 ±0.03 23 ±6 

0.77 ±0.03 20 ±5 

0.91 ±0.02 18 ±4 

1.07 ±0.02 11 ±2 

1.23 ± 0.03 4.3 ±1 

1.42 ±0.02 4.4 ±1 

1.73 ± 0.03 2.5 ±0.8 

1.94 ±0.04 2.2 ±0.8 

2.40 ±0.04 2.6 ±1 



0.58 ± 0.03 26 ± 7 

0.77 ±0.02 23 ±6 

0.90 ±0.02 25 ±6 

1.06 ±0.02 13 ±4 

1.21 ±0.04 4.4 ±2 

1.41 ±0.03 4.2 ±1 

1.75 ±0.05 1.8 ±1 

1.88 ±0.06 2.0 ±1 

2.36 ±0.06 2.6 ±1 



Additional Parameters 



r ••• 1.9f 

r ■ • • 3.0 ± 0.7 

X^/dof 25/26 46/49 



Note. — Errors are 90% confidence limits for 
one interesting parameter. Fixed parameters are 
marked with f. 

'^Energy of line center in keV. 

''Line flux in 10~^ photons cm~^s~^. 

'^Normalization of power law in units of 
10~^ photons keV~^cm~2s~^ at 1 keV. 
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Table 4 


. NGC 3393 Soft Spectra 


Nucleus 


Center 


^ line 


^ line 


^ Hne ^ line 


0.57 ±0.03 
0.70 ±0.03 
0.87 ±0.02 


26 ±7 
13 ±5 
12 ±3 


0.40 ± 0.09 91 ± 60 
0.56 ±0.01 98 ±20 
0.72 ±0.02 33 ±6 
0.88 ±0.01 33 ±5 


Additional Parameters 


kT (keV) 

AC 

thermal 

xVdof 


1.6 ±0.3 
3.3 ±0.7 
22/22 


1.7±0.2 
8.8 ± 1 
68/53 



Note. Errors are 90% confidence limits for one 

interesting parameter. 

''Energy of line center in keV. 

''Line flux in 10~^ photons cm~^ s~^. 

'^Normalization of thermal component in units of 
IQ-^K, where K = (IQ-^'^/i'iTTD^)) J UeUHdV, D 
is the distance to the source (cm), rie is the elec- 
tron density (cm~^), and nn is the hydrogen density 
(cm-3). 
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Table 5. Observed Flux 





-''0.5—2 


.f'2 — 10 




(10~-^^ergcm~^ s~^) 


(10~-^^ergcm~^ s~^) 


ESO 428-G14 nucleus 


0.91 ±0.2 


3.8 ±0.5 


ESO 428-G14 center 


1.4 ±0.1 


4.7 ±0.5 


ESO 428-G14 galaxy 


1.7 ±0.3 


6.5 ±0.7 


NGC 1386 nucleus 


1.1 ±0.2 


2.1 ±0.1 


NGC 1386 center 


1.6 ±0.2 


2.1 ±0.4 


NGC 1386 galaxy 


2.2 ±0.2 


3.8 ± 1 


NGC 3393 nucleus 


0.73 ±0.07 


2.1 ±0.04 


NGC 3393 center 


2.0 ±0.2 


2.4 ±0.07 


NGC 3393 galaxy 


2.7 ±0.2 


4.2 ±0.9 


NGC 5347 nucleus 


0.24 ±0.04 


2.2 ±0.4 


NGC 5347 northeast 


0.035 ± 0.01 




NGC 5347 galaxy 


0.34 ±0.1 


2.3 ±0.3 


NGC 7212 nucleus 


0.43 ±0.2 


3.9 ±0.9 


NGC 7212 center 


0.67 ±0.3 


5.4 ±0.6 


NGC 7212 galaxy 


0.92 ±0.3 


6.5 ±0.7 



Table 6. Observed and Intrinsic Luminosity 



Galaxy 


-^'063, 0.5-2" 






LAGN,Fe^ 


Lagn,pex'^ 


Lagn,[oiii]'^ 


Lir" 




r . f 


ESO 428-G14 


6.6 X 10=*'' 


2.3 X 10*" 


5.5 X 10^'' 


3.0 X 10*2 


3.6 X 10*1 


7.6 X 10*2 


4.1 X 10*3 


5.7 X 10*° 


1.3 X 10*3 


NGC 1068 


1.3 X lO"'! 


1.1 X 10"! 


1.6 X 10*0 


7.7 X 10*2 


6.3 X 10*1 


2.5 X 10*3 


1.0 X 10*5 






NGC 1386 


2.0 X 10-^'^ 


3.9 X lO^'' 


1.1 X lO^'' 


6.6 X 10*1 


5.2 X 10*" 


5.8 X 10*1 


1.3 X 10*3 


1.0 X 10*° 


2.6 X 10*2 


NGC 3393 


2.3 X 10*0 


7.5 X 10*" 


1.5 X 10*0 


7.2 X 10*2 


1.0 X 10*2 


1.7 X 10*3 


9.6 X 10*3 


1.5 X 10*1 


2.6 X 10*3 


NGC 5347 


3.2 X 1039 


3.0 X 10*° 


6.4 X 10^9 


3.0 X 10*2 


4.8 X 10*1 


6.5 X 10*2 


3.7 X 10*3 


9.5 X 10*° 


2.0 X 10*3 


NGC 7212 


9.1 X 10"° 


8.2 X 10*1 


1.2 X 10*1 


5.2 X 10*3 


1.3 X 10*3 


4.2 X 10*3 


5.4 X 10** 






Mrk 3 






2.0 X 10*1 


4.8 X 10*3 




5.5 X 10*3 


2.8 X 10** 







Note. — All luminosities axe in units of erg s~i. 

''Observed in nuclear spectrum, not corrected for absorption. 

''Intrinsic 2—10 keV luminosity, determined from Fe line luminosity. 

"^Lower limit on intrinsic 2—10 keV luminosity, determined from PEXRAV spectral model. 
"1 Intrinsic 2—10 keV luminosity, determined from [O in]A5007 luminosity. 

''8-1000/im bandpass, from IRAS flux. 
^Lower limit, from optical images. 
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Table 7. NGC 1068 Nucleus 



line 


line 


— u.uo 


350+^f 


0.55 ±0.01 


180+26 


o.7it°j; 


7ltl 






n QQ+0 03 
'-'•^^-0.04 


23_6 


1.12 ±0.02 


16 ±4 


1.27 ±0.30 


lOt^ 


1 no+0.07 

-'■•O»_o.i4 




1 t;o+0.18 

i-'J'J-o.io 


4.7 ±2 


-■-■00-0.04 


3.6 ±2 


1.88 ±0.02 


6.7 ±2 


2.20t°f 


3.3 ±2 


9 C1+0.09 


6.6t^ 


q no+5.7 
"^•"^"*-0.04 


2.7tl 



Additional Paramotcrs 



xVdof 45/42 



Note. — Errors are 
90% confidence limits for 
one interesting parame- 
ter. 

'^Energy of line center 
in keV. 

''Line flux in 
10~^ photons cm~2 s~^. 



Table 8. Emission Measures 



Galaxy (Region) 


EMovii 




elX 


EMHa 




(cm" 




(cm" 




(cm" 




ESQ 428-G14 (center) 


8.0 X 


1020 


2.6 X 


1021 


4.3 X 


1024 


NGC 1386 (center) 


1.2 X 


1021 


1.7 X 


1021 


2.1 X 


1023 


NGC 3393 (center) 


2.6 X 


1021 


3.0 X 


1021 


9.0 X 


1022 


NGC 5347 (nucleus) 


1.8 X 


1021 


1.7 X 


1021 


1.4 X 


1024 



-37- 



Table 9. Additional Sources 





R.A. 
(J2000) 




Dccl. 
(J2000) 


Count Rate 
(10-''cts-l) 


Distance 
(") (kpc) 


ESO 428-G14 
















7 


16 


25.22 


-29 


20 


35.19 


2.1 ± 1 


102 


12 


7 


16 


25.45 


-29 


18 


22.18 


1.8 ± 1 


100 


11 


7 


16 


26.41 


-29 


21 


3.43 


2.4 ± 1 


113 


13 


7 


16 


26.75 


-29 


18 


1.27 


2.1 ±1 


105 


12 


7 


16 


27.30 


-29 


18 


26.03 


2.0 ± 1 


81 


9.1 


7 


16 


97 fifi 


-29 


1 Q 


90 zlQ 
zu.^y 


1.1 ±1 




o.o 


7 


16 


29.23 


-29 


21 


20.45 


2.2 ±2 


114 


13 


7 


16 


29.95 


-29 


18 


20.32 


26 ±3 


70 


8.0 


7 


16 


30.07 


-29 


21 


22.19 


z.z ± i 


114 


13 


7 


16 


30.29 


9Q 


21 


17.81 


9 -1- 9 
o.z nz z 


110 


12 


7 


16 


30.63 


9Q 


18 


35.28 


o.o zc z 


54 


6.1 


7 


16 


30.80 


9Q 


19 


35.84 


o.o ^ z 


9 


1.0 


7 


16 


31.50 


9Q 


21 


12.66 


n -k 1 


104 


12 


7 


16 


31.85 


9Q 
-zy 


19 


52.35 


-\- 'i 

lO ^ o 


25 


2.8 


7 


16 


31.93 


9Q 
-zy 


19 


17.77 


A 9 -1- 9 


15 


1.7 


7 


16 


32.37 


9Q 
-zy 


17 


47.56 


12 rb 2 


103 


12 


7 


16 


33.01 


9Q 
-zy 


18 


1.59 


R 9 -1- 9 

U.Z ZL Z 


91 


10 


7 


16 


34.02 


9Q 
-zy 


20 


22.37 


9 0-1-9 

Z.U ZL Z 


65 


7.4 


7 


16 


37.06 


-29 


19 


28.. 36 


3 2 zt 2 


77 


8.7 


7 


16 


37.31 


9Q 
-zy 


20 


48.43 


"1 7 -1- "1 
-L. ( zn -L 


113 


13 


7 


16 


37.33 


9Q 
-zy 


18 


54.98 


OZ ^ 4: 


87 


9.9 


7 


16 


37.63 


-29 


20 


20.76 


2.6 it 1 


99 


11 


NGC 1386 
















3 


36 


37.36 


-oo 


59 


10.53 


LL III Z 


117 


7.0 


3 


36 


39.17 


-ou 


00 


46.09 


11 ill 2 


98 


5.9 


3 


36 


oy. uu 


-35 


oy 


OA 


4.9 it 2 


oo 


^ 1 

O.J- 


3 


36 


39.58 


-ou 


00 


56.63 


9 7 i 2 


100 


6.0 


3 


36 


43.56 


-ou 


00 


03.90 


1 fi + 

-LU ^ O 


32 


1.9 


3 


36 


44.18 


-36 


00 


16.92 


3.8 i 2 


31 


1.9 


3 


36 


46.31 


-35 


59 


48.56 


21 it 3 


9 


0.5 


3 


36 


46.50 


-36 


00 


05.04 


12 it 3 


9 


0.5 


3 


36 


46.75 


-35 


59 


46.56 


25 it 4 


13 


0.8 


3 


36 


47.32 


-36 


00 


24.99 


5.2 it 2 


31 


1.9 


3 


36 


47.61 


-35 


59 


00.08 


3.6 ± 1 


60 


3.6 


3 


36 


50.22 


-35 


59 


58.01 


19 it 3 


49 


2.9 


3 


36 


51.77 


-35 


59 


08.61 


7.6 it 2 


84 


5.0 


NGC 3393 
















10 


48 


15.67 


-25 


10 


0.24 


9.2 it 2 


107 


28 


10 


48 


17.66 


-25 


9 


37.24 


0.8 it 1 


79 


20 


10 


48 


19.13 


-25 


11 


27.01 


17 ±3 


119 


31 


10 


48 


26.01 


-25 


10 


23.36 


13 ±3 


53 


14 


10 


48 


26.03 


-25 


9 


15.24 


1.0 ± 1 


45 


12 


10 


48 


27.64 


-25 


8 


22.87 


16 it 3 


99 


26 
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Table 9 — Continued 



R.A. Decl. Count Rate Distance 

(J2000) (J2000) (lO-^cts-i) (") (kpc) 



10 


48 


30.61 


-25 


10 


6.05 


54 ±5 


100 


26 


NGC 5347 
















13 


53 


8.33 


33 


29 


35.45 


2.4 ± 1 


119 


19 


13 


53 


10.69 


33 


29 


22.41 


2.1 ± 1 


89 


14 


13 


53 


17.72 


33 


29 


8.25 


2.3 ± 1 


19 


3.1 


13 


53 


18.33 


33 


27 


38.84 


6.0 ±2 


109 


18 


13 


53 


23.85 


33 


30 


31.41 


6.2 ±2 


99 


16 


13 


53 


24.26 


33 


29 


36.85 


1.9 ± 1 


81 


13 


13 


53 


24.39 


33 


28 


38.76 


2.3 ±1 


96 


16 


13 


53 


26.38 


33 


30 


4.86 


6.7 ±2 


114 


18 


NGC 7212 
















22 


06 


54.36 


10 


14 


09.02 


9.1 ±2 


113 


63 


22 


06 


59.05 


10 


15 


40.47 


7.9 ± 2 


108 


60 


22 


06 


59.68 


10 


13 


13.03 


9.3 ± 2 


59 


33 


22 


07 


01.96 


10 


15 


16.16 


2.3 ± 1 


75 


42 


22 


07 


02.54 


10 


14 


16.08 


7.9 ±2 


17 


9.5 


22 


07 


04.23 


10 


13 


02.70 


11 ±2 


67 


37 


22 


07 


04.85 


10 


14 


22.53 


3.8 ± 1 


47 


26 


22 


07 


04.89 


10 


15 


08.99 


4.0 ±1 


80 


44 



Note. — Units of right ascension are hours, minutes, and seconds, and units 
of declination are degrees, arcminutes, and arcseconds. 
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Fig. 1. — X-ray image of ESO 428-G14. The image has been adaptively smoothed and is 
scaled hnearly. The contours are spaced logarithmically. 
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Fig. 2. — Continuum-subtracted Ha image of ESO 428-G14 with linearly-scaled X-ray con- 
tours overlaid. The similar morphology of the X-ray and optical hne emission suggests their 
common photoionized origin in all the sample members. 
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Fig. 3. — The prominent Fe Ka lines are evident in all the observed nuclear spectra. The 
best-fitting model of the hard emission, which includes a pure reflection continuum and the 
emission line, is plotted with the data. We show binned data here but use the unbinned data 
in the model fitting. 
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Fig. 4. — Nuclear spectrum of ESO 428-G14 {crosses) and best-fitting model {red histogram), 
{upper panel) Observed spectrum and total model, uncorrected for detector sensitivity, {lower 
panel) Emitted spectrum and model, corrected for detector sensitivity. Individual model 
components are plotted with dot-dashed lines. The hard X-ray spectrum shows the reflected 
AGN continuum and the Fe Ka line, while 9 unresolved Gaussian line profiles account for 
the soft X-rays. 
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Fig. 5. — Ratio of intrinsic AGN 2-10 keV luminosity determined from the Fe Ka line 
luminosity to intrinsic AGN 2-10 keV luminosity determined from the [O III] luminosity vs. 
[O III] luminosity. The luminosity is in units of erg s^^. The sample members are plotted 
as crosses, except NGC 1068, which is plotted as a triangle. Two additional Compton thick 
galaxies from the literature that do contain circumnuclear starbursts (NGC 5135 and NGC 
7130) are plotted as diamonds. The line is not a fit to the data and shows the good agreement 
of the two methods for determining intrinsic luminosity. 
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Fig. 6. — Nuclear spectrum of ESO 428-G14 fit with a thermal model {kT = 0.7 keV) in 
addition to the hard reflection components. The model flts the soft X-ray emission poorly, 
producing a broad emission hump around 0.9 keV, rather than the observed narrower fea- 
tures. 
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Fig. 7. — This ACIS spectrum of the nucleus of NGC 1068 shows narrow emission features, 
similar to the the ACIS spectra of the nuclei of the other galaxies. High spectral resolu- 
tion observations of NGC 1068 confirm the dominant role of photoionization in producing 
the observed emission. The observed spectrum [crosses) and total model [red histogram), 
uncorrected for detector sensitivity, are plotted. 
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Fig. 8. — Observed soft (0.5-2 keV) vs. hard (2-10 keV) luminosity, uncorrected for absorp- 
tion, in units of erg s~^. The hne is the relationship for a completely unobscured AGN, not 
a fit to the data. These Compton thick AGNs would appear to be unobscured and weak 
based on their hardness ratios alone. Symbols as in Figure 5. 
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Fig. 9. — X-ray "color-color" plot, based on ratios of counts in the soft (0.3-2 keV), medium 
(2-5 keV) and hard (5-8 keV) bands that Chandra would observe. Each curve repre- 
sents a different AGN continuum model. Solid curves show fixed AGN obscuration with 
Nh = 3x 10^^ 102^ 3 X 102^ and lO^^cm'^ (from top to bottom). Ahernatively, the AGN 
contribution is calculated from the pure reflection model, either excluding {dotted line) or 
including [dashed line) an Fe Ka line having EW = 1 keV. An empirical photoionization 
model is mixed with each of these AGN components. Along each curve, crosses are plotted 
where the photoionization model contributes 0, 20, 40, 60, and 80% of the total observed 
counts (from left to right). Thus, only the leftmost point of each curve describes a pure 
AGN model, one without any reflected or scattered soft X-ray component. A constraint on 
the efficiency with which the AGN continuum is converted into line emission excludes higher 
values of the soft/medium ratio in the moderate absorption models, distinguishing them 
from the reflection models. Only Compton thick AGNs are found in the region to the right of 
the thick solid line. Such colors would require excessively high reprocessing efficiency when 
absorption is Compton thin, with Lphotofi.5-2/ Lagn,2-io > 0.03. These strongly excluded 
regions of the moderately absorbed AGN models are also marked in red. Green identifies 
mixes where the reprocessing efficiency is high, with 0.01 < Lphoto,o.5-2/ Lagn,2-w < 0.03, al- 
though not strictly requiring Compton thick obscuration. The colors of galaxies of this study, 
the comparison Compton thick AGNs that also contain starbursts (both filled symbols), and 
several others that are not Compton thick [open symbols) are plotted. 
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Fig. 10. — Two scenarios to observe the intrinsic AGN continuum {solid lines) in reflection, 
where an obscuring torus (gray) blocks direct views. For an observer at the right of the 
cartoon, the reflecting medium may be the far side of the torus, and the reflected hght is 
not further obscured {dotted line). Ahcrnativcly. the blocking medium may also serve as 
the reflecting medium, in which case the reflected spectrum itself is obscured {dashed line). 
The unobscured continuum reaches the material of the photoionized emission region {pale 
gray). This figure is not to scale. The torus inner radius is on the order of 1 pc (the dust 
sublimation radius), while the photoionized emission extends over hundreds of pc. 
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Fig. 11. — Ratio of total IR luminosity to intrinsic AGN 2-10 keV luminosity determined 
from the Fe Ka line vs. intrinsic AGN luminosity determined from the Fe Ka line. The 
luminosity is in units of erg s~^. These Compton thick AGNs that do not contain starbursts 
(crosses) show IR/Lagn ratios similar to unobscured AGNs (dotted line). The galaxies 
that do contain strong nuclear starbursts (diamonds) and NGC 1068 (triangle), which has a 
kpc-scale starburst, exhibit significant excess IR emission. 
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Fig. 12. — Ratio of mid-IR luminosity to intrinsic AGN 2-10 keV luminosity determined 
from the Fe Ka line vs. intrinsic AGN luminosity determined from the Fe Ka line. The 
luminosity is in units of erg s^^. Similar to Figure 11, the luminosity ratios of the Compton 
thick AGNs that do not contain starbursts {crosses) are similar to those of unobscured AGNs 
{dotted line). While the mid-IR excesses of the starburst galaxies are less severe than their 
total IR excesses, mid-IR data alone do not accurately measure AGN luminosity when stellar 
contamination is significant. 



